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ABSTRACT 


The analysis of the solutions of large scale optimize- 
tion models is very difficult without computer aids because 
typical outputs may exceed 40,000 lines. A computer system 
to allow efficient storage and interactive analysis of the 
solution file was designed and implemented in FORTRAN. 

This system, called ATHENA, has been initially designed 

for problems with up to 30,000 rows and columns; tests show 
the system has fast response time for these large problems. 
A description of the system, its data structures and 


commands as well as a user's manual is included. 
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i INPRODUCTION 


The generation, solution and analysis of large scale 
mathematical programming models presents Significant problems 
in data handling. For a typical large scale model the output 
from the computer may exceed 40,000 lines and it is very 
difficult and time consuming for the user to go through that 
much paper to extract the information he needs. Space is 
also needed to store all these outputs and it is difficult 
to provide routine access to old solution files. 

Over the past several years, faculty and students at the 
Naval Postgraduate School and the University of California 
at Los Angeles have been cooperatively developing theory and 
algorithms to solve large scale linear, nonlinear and 
integer optimization problems. This research has been 
furthered by the development of several software systems to 
solve large scale optimization problems. ATHENA is part of 
that development and was built to satisfy a pressing need 
to be able to quickly and easily analyze solutions to these 
large problems. 

The research in large scale optimization at the Naval 
Postgraduate School has concentrated on providing economic 
solutions to current Department of Defense problems. One 
such project that was going on concurrent with the development 
of ATHENA was a large, medium range capital budgeting problem 


mmiat required the solution of a mixed integer programming 





problem with 11,687 constraints and variables [6]. ATHENA 
was used successfully with this project, and performance of 
ATHENA on this problem 1s reported below. 

System ATENA has been developed to handle output from 
these large optimization problems by enabling the user to 
get the information he needs interactively through a com- 
puter terminal and by economically storing the large files 
in packed form in low cost media. The features of the 
system may be summarized as follows: 

1. Quick and accurate answers to simple guestions that 


are tedious and error prone to address manually, e.g., 


How many of a set of variables are = 1.0? 
Or How many are greater than Q? 
or What variables are in a specified range? 
or What constraints are satisfied exactly? Etc. 
2. In large scale mathematical programming models the 


names of rows and columns are constructed systematically so 
that groups of variables with relationships in the real 
world have similar names. Using the system one can have 
automotic, easy and accurate answers for many interesting 
properties of these groups. (For instance, the average 
value of all the variables whose names begin with X, etc.) 

3. ATHENA can also be used as a basis for a simple 
report writer. 

4. The system provides very compact computer storage: 
maewOonrloqmmal solution file is typically packed into 1/10 


of its normal volume. For example, a solution file from the 





IBM MPS/360 package of a linear programming model with 
12,000 rows and columns occupies 1.5 magnetic tapes 2400 
feet long at 800 BPI in original unblocked form, but in 
packed form the solution occupies approximately 31 feet of 
magnetic tape. 

9. There 1s systematic and economic file organization 
with easy and accurate access. 

See etl lewcemucture for multiple runsMof the problem 
for comparisons is available. 

7. The system requires modest resources (core, time) 
in a time sharing environment for a large amount of data. 

8. The system has portability and allows easy change 
or expansion since the whole system has been developed 
in FORTRAN. 

ATHENA waS inspired by a similar system developed for 
the Department of Energy by O'Neil and Sanders [5] called 
PERUSE. PERUSE was developed to aid in the analysis of large 
linear programming energy models. A study of the needs of 
Naval Postgraduate School students and faculty showed that 
additional capacities beyond those in PERUSE were necessary 
to support current and future research in large scale 
Seeimezation. In addition to the standard MPS output, it 
was determined that ATHENA should support the experimental 
weimizadeilon system XS [3]; the output of this system con- 
tains in addition to standard MPS output, upper and lower 
penalties that implement the ‘elastic formulation' of linear 


models that is unique to XS. ATHENA also had to support the 


10 





use of a preprocessor PREP [2] that reformulates the 
Original optimization problem to an equivalent reduced 
problem with fewer rows and/or columns. A study of past 
and current modeling efforts at the Naval Postgraduate 
School identified additional commands that would help in 
the analysis of large models. 

ATHENA was designed to be as much as possible a direct 
extension of PERUSE. Almost all the commands and options 
of PERUSE have been included with the identical names and 
syntax whenever possible. A summary of the extensions is 
listed under the section LIMITATIONS - EXTENSIONS, of the 


moer Ss Manual. 


Jib 








Eis ~SoOriyere DESCRIPTION 


A. GENERAL 
The whole system consists of 3 basic subsystems (see 
Bagure 1). 
t. PARSER - Code Generator 
This subsystem accepts aS input a Query, parses it 
examining the syntax according to the productions of the 
Query Language and generates the corresponding internal 
code or gives information for syntax errors. The internal 
codes for each Query are shown in the program list. 
Ae SBR Eg) Dios 
This subsystem accepts as input either (1) an 
unpacked solution file in 'standard' format which it packs 
and saves for future reference, or (2) an L.P. solution 
file in packed form from a previous session. 
3 Iles Has 
This subsystem, uSing the code generated from the 
PARSER, searches the packed solution file and prints out 


the information requested. 


B. LANGUAGE 
IPeenost) sanguage 
The system has been developed in a portable subset 
of FORTRAN IV. FORTRAN was chosen for the following reasons: 
a. FORTRAN is a general language available at almost 


any computer installation, so the system can be used with 


ps 





~~ 


QUERIES 
_ PACKED 
PARSER 
CODE INTERPRETER | 
GENERATOR 


FIGURE 1. System ATHENA 


13 





any contemporary hardware. Non-IBM systems may require 
some program modifications. 

b. Since FORTRAN is a high level language, the 
development time was low. 

c. The response time for each Query is acceptable 
and there is little need for faster responses or enhanced 
Pee rcLency . 

d. Extensions and changes of the system can be 
eaSily implemented. 

e. There is good system support for FORTRAN. In 
particular, ATHENA was developed with the FORTRAN G compiler 
oo LBM 360/67. 

2 QUE Language 

The set of acceptable Queries is divided into two 
main categories: 

a. Control Queries 

Control Queries provide commands to the system 
to perform specific tasks, but usually do not use the 
sOlution file. Examples of Control Queries are those that 
accept comments for self documentation of the output, 
print headings for the output, terminate the use of system, 
Se ¢ 

b. Command Queries 

Command Queries use the solution file to extract 
the information asked for. Each Command Query consists of 
various fields separated by at least one blank. Some of the 
fields are optional, while others are required. An internal 


code number is generated by parsing the Command Query for 





each field depending upon the analysis of that field and 
the previous fields in the Command Query. 

The code generated by the PARSER is executed 
by the INTERPRETER which consists of a set of progams 
(subroutines) activated by the code numbers, to get the 


required information from the file. 


ee DATA STRUCTURES 
1. Preview 

There are some observations about the solution file 
of a linear program, especially of a large scale one, that 
lead to the use of a special data structure for storing a 
solution file in less memory space than it would otherwise 
require. 

For each row or column the following information 1s 
usually included in the solution file. 

NUMBER of row or column. 

NAME, usually 6-8 alphanumeric characters. 

Sramesemusually 2 characters, e.g., BS for BASIC, 
meeeror LOWER LIMIT, etc. 

AGM OLY LEVEL, Lor each row or column. 

SieeineacTiIViTY for rows or INPUT COST for columns. 

LOWER LIMIT 

UPPER LIMIT 

Denes ert vil. for rowS or REDUCED COST for columns. 

UPPER PENALTY and 

LOWER PENALTY for the elastic linear programming 


Syecem, XS [3]. 


als: 





There 1s redundant information in each record. Some 
of these redundancies are the following: 

The explicit number for each row or column may be 
represented implicitly by the ordinal position of the row or 
column in the file. Rows usually precede columns in solution 
files. 

A large number of the ACTIVITY LEVEL values will be 
zero. The same is true for the SLACK ACTIVITY, LOWER LIMIT, 
DUAL ACTIVITY and PENALTY values. 

In many cases there will not be LOWER or UPPER LIMITS 
S@ee i NALTIES. 

PENALTIES in some cases may be infinite. 

Hihemethe Status Of a row or column is ‘fixed’, then 
ACTIVITY LEVEL, LOWER and UPPER LIMITS are all the same 
number. 

Each row or column can be in only one of its 
possible states. 

Moreover, analysts who have experience with large 
scale Linear Programming have observed that most of the 
numbers of the solution file are the same. For example, 
most of the numbers for LIMITS are the same for a large 
number of rows or columns. For purposes of analysis, it 
is rarely necessary to have more than five decimal digits 
of precision for problem values. Indeed, some large prob- 
lems cannot be solved with even this degree of significance. 


Accordingly, IBM single precision REAL*4 representation 


INS 





is adequate for our purposes. Conversion to REAL*8 extended 
precision requires trivial program modifications. 

Based on the above observations, two types of data 
Structures for storing the solution file have been developed. 
The first one (SPARSE) exploits the redundant information 
into each individual record. The second (SUPERSPARSE) takes 
advantage of the last observation by storing each distinct 
number only once for the entire file. It is the responsibility 
of the user to select the data structure type that 1S appro- 
priate for each solution file. SUPERSPARSE is probably 
Superior for problems with less than half of all numbers 
possessing distinct real values. What follows is a detailed 
description of these two data structures. 

Zo) SEARS® Data Structure 

PMacmenterse SOluttlon file 1s stored in contiguous 
memory (8-bit bytes) as a one-dimension array called SOLFIL, 
in the following way: 

a. The first 16 bytes (Four 4-byte words) are used 
to keep information for: 

(1) The size of the file in 4-byte words. 

(2) The type of data structure used to pack the 
file (SPARSE or SUPERSPARSE). 

(3) The number of rows and columns of the file. 

b. For each row or column, 12 sequential bytes are 
required, organized as follows (see Figure 2}. 

(1) The first 8 bytes hold the name of the row 


@empecolumn left justified, one character per byte. 


Baer 





S—-BY TE Wali i 


225.15 BIT MAP 
| 2-BYTE POINTER =BwiE; POINTER 


— 8-BYTE NAME-2 





2-Svtee oe. 


eet POINTER | BYTE POINTER 
8S-BYTE NAME-n 
: WAP 
ie 2-BYTE POINTER 


4-BYTE VALUE 










4-BYTE VALUE 


4-BYTE VALUE 


f—-EB i TE WALUE 





4-BYTE VALUE 


| 4-BYTE VALUE Ete VAGUE 


Seemed Sige ye 


ea — sais ay Lor veo 
Aen wer y AU 
f—pr le VAnuUe 


FIGURE 2. SPARSE Data Structure 
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(2) The next 2 bytes are used (as 16 bits) to 
represent various characteristics associated with that row 
eaecolumn. 

(3) The last 2 bytes are used aS a pointer to 
the first number stored from the current record. 

c. The 16 bits from (2) above are organized in 4 
groups of 4, 7, 4 and 1 bits, respectively, taken from higher 
mowower order. 

The first group of 4 bits represents the status 


@e the current row or column. 


syligk PATTERN STATUS 

C200 0 EN SPP EAS iE Gs ) 
CoO yi. BS (Basi) 

OF Oi.L 0 LL (LOWER LIMIT) 
Ceor lt 1 Uwe E ER bitin) 
Cas 0 0 EO Gs RED) 


The following status indicators are reserved for use with 


the program PREP [2] 


ori 0 1 VG VOrD -CeremnN) 

CE a @, SC (SINGLETON COLUMN) 

C) 1c eae ag FC (FIX COLUMN) 

io 0? 0 BC (BOUND CHANGED) 

aL OmaO ced VR (VOID ROW) 

EO 1 0 SR (SINGLETON ROW) 

iO 1 1 RR (REDUNDANT ROW) 

mt oO 0 FR (ROW FIXES VAR. AT BOUND) 
es OL ER (DOUBLETON EQUATION) 


MS, 





Jey Jae 6 


me. 1 


TR (TIGHTEN RANGE) 


PP Reserved for PREP [2] 


The next group of 7 bits represents the characteristic of 


zero or nonzero values for the record. 


oe NO: 


iL 


10 


0 


BIT VALUE: 


CHARS CLE RISTIC 


2G Ee seEVEL NONZERO 


AGLiViEyY LEVEL ZERO 


SLACK, COST 
SLACK GOST 
OWE he bleak i 
LOWER LIMIT 
UPPERS TMT 


ONS ods, EVO sG ah 


NONZERO 


ZERO 


NONZERO 


ZERO 


NONZERO 


ZERO 


DUAL/RED. COST NONZERO 


DUAL/RED. COST 


LOWER LIMIT 


LOWER LIMIT 


CPP ER EaaT 


UPPER saeaiT 


ZERO 

0, Go Me: 
DOESN: T ESD 
EAISLS 


DOESN. T«bxXtST 


The next group of 4 bits represents the characteristics for 


PENALTIES. 


See PATTERN 


00 0 0 


Orr O a a. 


0010 


FO? lek 


UPPER ONLY 


i StS. 


ZERO 


ZERO 


NUMBER 


Z 


LOWER PENALTY 
ZERO 
NUMBER 
eve LN Ey 


CERO 


0 





oreo 0 NUMBER NUMBER 


Clk 0 )aa NUMBER LNESEN TTY 
Oro NE ZERO 

Sg ae Eo UNEINITY NUMBER 
iO. 410 0 Chee Ty Then. Py 


The rest of the bit permutations are not used. 

The last (0 bit) is used by the interpreter to 
mark the active and nonactive records when the user uses 
the ACTIVE or DEACTIVE commands to avoid searching of the 
entire file. 

All the above groups of bits are stored together 
as a 16 bit binary number, which is stored in 2~byte half- 
word. ATHENA has provisions for the use of 16 bit halfwords 
representing absolute magnitudes of 0~ 65535, and can extract 
any component bits of the halfwords as necessary. (In this 
sense, the usual signed magnitude of IBM/360 halfword 
integers 1s ignored.) 

d. The (nonzero, noninfinite) number values which 
must be stored are located immediately after all the informa- 
tion above. If the file represents a problem with M rows 
and N columns, then location INDEX - where INDEX=(N+M) *3+4+1 - 
of the SOLFIL array is the first eligible location for 
storing number values. The value of INDEX is kept ina 2- 
Peer pOinter associated with each row and indicates for that 
row the location of the first value stored. The sequence for 


storing these numbers for each row is: 


Zi 





ACTIVITY LEVEL, SLACK/INPUT COST, LOWER LIMIT, UPPER LIMIT, 
DUAL/REDUCED COST, UPPER PENALTY ,~ LOWER PENAMGTY . 
o-) SUPERSeaRes Data Structure 

This type of data structure takes advantage of the 
fact that in most problems many number values in the solution 
file are the same. Each distinct value is stored only once 
and a 2-byte pointer is used to access this value when 
needed. This is the only difference from the SPARSE repre- 
sentation (see Figure 3). 

The array with the packed solution file is now 
separated into 3 parts: 

a. The first part is exactly the same as in SPARSE. 

b. The second part is substantially the same with 
the following differences: 

(1) It consists of 2-byte halfwords instead of 
4~byte words. 

(2) Each halfword is a pointer to the third 
part of the array where the distinct number values are 
stored. 

Ceeie caine Pateeeensises OL a pool of 4-byte 
words, each representing a distinct real number value. The 
pointers to the distinct real number values are relative 
addresses in the real number pool, so. a file which 1s packed 
with a different array size can be used with the current 
pointers providing the array size is large enough to holc 


the file. 


Ze 





— i 





Ceey fees | 


Zee eis LL MAP 


25s POINTER 
Se bY TE NAME-2 
RRS EE 2 al =proe Bit MAP 
— =—2 ue | 2=BYTE POINTER  _—s 
S=by iE ei NAME=n 
| « 2=BYTE BIT MAP ss BYE | « 2=BYTE BIT MAP ss MAP 
| -« 2=BYTE POINTER eke Eee POT Niet 


Sea =e | ss 2=BYTE INDIRECT PTR. | eS 
ie Ere ser. oes, INDERECT PER. 


— ee Re Sole eID iRser PTR. 





Zoe | ss 2=BYTE INDIRECT PTR. | ER. 
Fe pills | 2-BYTE INDIRECT PTR. | POR: 


ee — ee otoriNer VALUE 
fp Lae Diol ENG VALUE 
4—-BYEEeDESTINGs VALUE 


2B ys DISTEINGY VALUE 
(HEY tee plot ive l VALUE, 


FIGURE 3. SUPERSPARSE Data Structure 
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Tift. USER'S MANUAL 


ae LNTRODUCTION 

The system ATHENA is a set of programs which accepts as 
input a linear programming solution file, packs it ina 
special data structure and interactively extracts specific 
information from that file through a set of Queries. 

The size of memeory which is required to run the system 
depends on the size of the file to be accommodated, and thus 
on the size of the original optimization problem. The user 
extracts information from the solution file with a Query 
Language, asking questions related to the solution of the 
problem represented in the file. 

The entire system has been developed in FORTRAN language 
for portability and better coordination with other Linear 
Programming procedures which are also written in FORTRAN. 

The Queries are self-documenting and their syntax follows 
closely the syntax of the English language. To avoid typing 
effort for experienced ATHENA users, short forms of Queries 
are provided. Only the characters comprising the short 
forms are interpreted by the system, with all subsequent 


contiguous nonblank characters ignored. 


B. QUERY LANGUAGE 


The Query Language consists of three subsets of Queries: 
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1. The Control Queries: 

With this subset of Queries the user controls 
mainly the output of the system, inserting comments, 
headings, etc. 

2. The Command Queries; 

With these the user communicates with the solution 

file and extracts the specific information he needs. 
5. Theyse! Command: 

This command qualifies the ATHENA queries to access 

only a subset of the problem file. 
1. Control Queries 

ay) Very SeSsnOres Porm. V) 

All the following Queries will be displayed with 
the output. This Control Query is useful when the OFFLINE 
printer 1S used for the output instead of the terminal, 
or when the system is used under Batch Processing; in these 
cases answers are transmitted to the output device without 
the corresponding questions if the system is not in VERIFY 
mode. 

See NOVER Thy (NOV > 

The following Queries do not appear with the 
output. This Control Query is most frequently used when a 
terminal is used for all output. The DEFAULT mode of the 
ATHENA system 1S NOVERIFY. 

c. PROMPT (P) 


The system responds with the prompt: 


' INPUT A COMMAND ' 


aS 





whenever it is ready to accept a Query. PROMPT is a 
DEFAULT mode of the system. 
d. NOPROMPT (NOP) 

Used to avoid the prompting phrase in the 
output, especially when the OFFLINE printer or the Batch 
Processing is used. 

e. H Any Character String 

When the first column of a Query is the letter 
H, then the character string is printed as is in the output. 
H 1s used to insert comments or headings in the output. 

any Character String 

When the first column of a Query is the character 
*, no action takes place. This is considered as a comment 
and 1S ignored. * is useful to insert comments and/or 
headings on the terminal output. 

Cpe eee) 

Used to end the current session. 
2. Command Queries 

A Command Query consists of several fields. Some 
fields are required and must always appear in a Command 
Query and others are optional. Each field is separated 
from the others by at least one blank character. The number 
of blanks between fields is not significant and a Query may 
Start at any character position in the command. The length 
of a Query cannot exceed 80 characters including the spaces 


between the fields. 
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The possible fields that can be included in a 


Command Query are: 


a. TYEE 
low So uLeea 
c. MASK 


do COND? TON 
See Or T LON 
The fields in a Command Query must appear in the 
above sequence and the first 3 of them must always appear, 
With only 2 exceptions. A detailed description follows of 
each individual field and the way that it may be used. 
ae Dee fel d 

This is the first field of the Query and may 
Start at any character position. This field can be one of 
the following: 

GomeDISPEAY {D) 

Used when all the records which meet the 
requirements of the other fields are to be displayed in the 
output in the sequence they are encountered starting from 
the beginning of the solution file. 

What portion of each individual record 
will be displayed depends on the PRINT OPTION field. 

Cee COUNT a(.GC,) 

Used when only the number of records which 
meet the requirements of the other fields is desired. COUNT 
is especially useful immediately preceding a DISPLAY command 


so the user will know in advance the size of output, avoiding 


2g 














unpredictably extensive printouts. For this Query the 
PRINT OPTION is ignored as meaningless. 
(3) ADD (A) 

Used when some numerical quantities of 
the qualified records are to be summed. The names of the 
numeric quantities of each record that will be added are 
given in the PRINT OPTION field. If no PRINT OPTION 
appears, all the numberic quantities of each record are 
added and their sums are displayed with appropriate labels. 

Since it is mathematically meaningless 
meracdd LOWER or UPPER LIMITS, or PENALTIES, they can not 
be summed or displayed. 

(4) AVERAGE (AV) 

Used exactly as the ADD command to display 
arithmetic averages. The sums are divided by the total 
number of the qualified records. 

(mee rlivATE (AC, (Syn. ACTIVE) 

With the ACTIVATE command the user can 
indicate a subset of the records of the solution file with 
Specific qualifications determined by the other fields so 
that subsequent Queries will implicitly refer only to that 
Subset. The user can expand the initial subset by using 
the ACTIVE command repeatedly to add new records to the 
active subset. 

The command ACTIVATE can minimize the 
Searching time for the required information in the active 


subset. Each time the ACTIVATE command 1s issued, the system 
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responds with the number of records added to the active 
Subset and the current total number of active records. 
(6) DEACTIVATE (DE) (Syn. DEACTIVE) 

Used to delete records with specific 
qualifications from the current active subset - created 
by the ACTIVE commands - or to eliminate any active file. 
The system responds with the number of records deactivated 
and the total number of records cceeree active. 

The entire active file can be deactivated 


Dy : 
_ DEACTIVE ALL or DE A‘. 


With this Query all the currently active records will be 


deactivated and the message: 
Pewee Fr Awe TL VE * 


will be printed out. Subsequent queries will refer to the 
entire solution file. 
Beewoemner field 

This field is mandatory and specifies whether 
the qualified records are ROWS, COLUMNS or BOTH. It may 
consist of one of the following: 

(1) ALL (A) 

Specifies that the entire file must be 

searched for the qualified records starting from the first 


ROW and continuing to the last COLUMN. 
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(2) COLUMNS (C) 

Specifies that the COLUMNS only will be 
searched for the qualified records starting with the first 
COLUMN and continuing to the last COLUMN. 

(a yreeeROWs = (RR) 

Specifies that the ROWS only will be 
searched for the qualified records starting with the first 
ROW and continuing to the last ROW. 

c. MASK field 
Specifies that any record is qualified for 
processing if the name of the record fits the MASK field. 
The MASK field is left justified and may contain 1 to 8 
characters. All the right unfilled positions up to 8 
characters are assumed to be the character *. The MASK is 
matched against the name, Starting from the left, character 
by character. Any character in the name is matched with 
a * in the MASK field. The MASK field is mandatory. 
EXAMPLES 
WeeiroemaAGk ) |’ X*****x*y' specifies all the names 
Starting with the letter X and having as the 8th (last) 
character the-letter Y. 
ie THe MASK "K'srs equuvalentewith the MASK 'X*******" 
and means all the names starting with the letter X. 
iii. The MASK '******yY' specifies all the names 
ending with the letter Y and it is NOT equivalent with the 


MASK "Y'. 
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iv. The MASK '‘'ABCDEFXY' specifies only this 
name and since the names of ROWS and COLUMNS are assumed 
to be inclusively smdenr = the searching of the file stops 
when the first match is made. 

v. The MASK '*' specifies ALL the names and 
may be used when no particular mask is desired. 
ad. CONDITION field 


The syntax of this field is: 
FOR ( <conditional phrase> ) 


The word FOR, left parenthesis and right parenthesis must 
always appear when the CONDITION field appears in a Query. 
There are two kinds of conditional phrases: 
The simple conditional phrase and the compound conditional 
phrase. 
(1) Simple Conditional Phrase 
There are 3 kinds of simple conditional 
phrases: The Relational, the Status and the Bound simple 
conditional phrases. 
(a) Relational Simple Conditional Phrase 


The syntax 1s: 
<Argl> <Relop> <Arg2> 


where Argl, Arg2 and Relop are one of the following: 


ies A Gl 
x Gris PGA Vvaeer lis el, 
Sor C IONE SEG ACTIVITY Or MIPUT COST 


Syl 








L icons LOWER LIMIT 


U Ow UPPER EMirt 
D ens DUAL eSACEiawryY Or REDUCED COST 
P Otc UPPER PENALTY 
W Bove LOWER PENALTY 
11. Relop 


Relational operators EQ, NE, GT, GE, LT, LE with the same 
meaning aS in FORTRAN. (Note, however, that there are not 
imbedding decimal characters as in FORTRAN.) 

ianibes cAaee2 

Arg2 is defined exactly as Aral 
with the enhancement that Arg2 may also be any integer or 
real number. Aro2 cannot be expressed as a floating point 
number in exponential notation. 
(b) Status Simple Conditional Phrase 


The syntax is: 


Siac bag Or oF <Flag- 


where Flag is one of the following: 


BS 1G BASLC 

LL EOws LOWER LIMT 

UL £Ox, UPPER LIMIT 

EQ 1g FIXED 

VC maha VOID COLUMN 

Sc Ou SINGLETON COLUMN 
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ne iene FIXED COLUMN 


BC Ox BOUND CHANGED 

VR Oe VOID ROW 

SR ews SINGLETON ROW 

RR fer REDUNDANT ROW 

ER For FREE ROW 

ER Eoss DOUBLETON EQUATION 
Ab Re a)1g TIGHTEN RANGE 


(c) Bound Simple Conditional Phrase 


The syntax is: 


<Argl> MINIMUM or <Argl> MAXIMUM 


where Argl is specified as in Relational Simple Conditional 
Phrase. The words MAXIMUM or MINIMUM can be abbreviated as 
MAX or MIN, respectively. This is used to extract those 
records which have the MAXIMUM or MINIMUM value in the speci- 
fied field with the specified MASK. The system responds 
with the first record encountered with the maximum or minimum 
value associated with Argl, and the total number of records 
that meet the requirements. This phrase may not be used 
Polen Aclivye Or DEACTIVE options in the TYPE field of the 
Query. 

(2) Compound Conditional Phrase 


The syntax of this phrase is: 


<Relational Cond. Phrase> <Log. Oper.> <Relational Conc. 


Phrase> 


Bie. 








or 
<Relational Cond. Phrase Log. Oper. Status Cond. Phrase> 


where Log. Oper. 1s OR or AND with the meaning of the corres- 
ponding logical operators. Note that the Bound Simple 
Conditional Phrase is not compatible for use in a Compound 
Conditional Phrase, since it exhibits no boolean value. 
Also, the Status Conditional Phrase must always appear after 
the logical operator. 

The CONDITION field as a field must be 
separated by at least one blank from the other fields of the 
Query. The word FOR, the left parenthesis, and the first 
element of the conditional phrase do not require separation 
by blank characters, nor do the last element of the condi- 
tional phrase and the right parenthesis. 

The CONDITION field is optional and need 
not appear in the Query. If it is not present, any record 
is qualified if the MASK field is satisfied. Using the 
ACTIVATE and DEACTIVATE commands the user can actually have 
unlimited length conditional phrases, by adding qualified 
subsets of records in the ACTIVE file. 

e. PRINT OPTION field 
Mais is the last field of a Query. It 1s 
optional, and if it does not appear the entire record which 
Satisfy both the MASK and the CONDITION fields are printed 


Out. 


34 





The elements of the PRINT OPTION field may be 


any combination of the following: 


1k ce OF Se, iby teh, BLY P, W 


with meanings as described in the CONDITION field. The 
Sutput will include information described in the PRINT 
OPTION with corresponding headings. The elements of the 
field can be separated by any number of blanks, by commas, 
SieemoOteyat ali. Tf both C and S appear in the PRINT OPTION 
neither of them is printed out. For the commands ADD and 
AVERAGE the default PRINT OPTION is X, C or S, D Since there 
is no meaning for LIMITS and PENALTIES. 

For all Queries that potentially require more 
than one output record for the answer (1.e., all Queries 
ereept COUNT ACTIVE, DEACTIVE and SET), the output will 
include the following entries for each record: 

NUMBER, NAME, STATUS and the entries specified in the PRINT 
OPTION field in the sequence in which they appear. At the 
end of the answer output for each Query the total number 
of qualified records is given. The heading for the output 
is determined by the SELECT field. If for this field the 
Option ALL is used, the heading will be the one for ROWS 
although COLUMNS may also be included in the output. 
See new on. Command 
By default each time a Command Query 1S issued the 


whole solution file is searched starting at the first ROW or 
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COLUMN and continuing by examining sequentially all the 
Becords. 

Queries may sometimes apply only to a small part 
of the solution file or to records whose relative position 
in the file is known. In these cases the SET command can 
cause searching to be initiated at a particular entry in 
the file and continued to another particular entry. Also 
a fixed step size can be specified for the search. Thus 
much computational effort can be avoided. 


The syntax for the SET Command is: 


SET <number l> <number 2> <number 3> 


where; 
number 1 is the number of the starting record 
number 2 is the number of the record to stop 
searching 


number 3 is the step for searching. 


All these numbers must be integers separated by at least one 
blank and the presence of all of them is required. These 
numbers alSo must be in the range of total number of recorecs 
por eae file. The SET limits apply to qualify any subsequent 
search of the file even if ROW or COLUMN subsets are speci- 
fied by a Query. 

EXAMPLE 

Suppose the solution file has 300 rows and 2000 


columns and the following SET command is issued: 
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See eno Loc0ny LO 


For all subsequent Queries: 

If the SELECT field of the Query is ALL then the searching 
Starts at the 18th record and continues through the 1500th 
record with step 10 (i.e., Record numbers 18, 28, 38, 
are examined). 

If the SELECT field 1s ROWS then the searching will start 
at the 18th row through the last row (300th) with step 10. 

If the SELECT field is COLUMNS then the searching will 
mearceat L8th column through the 1500th column (or equiva- 
Menmely the 318th record through the 1800th record) since the 
number of columns is greater than 1500. 


To restore default settings, use: 


oh Lari Cite. ok masD. 


C. ERROR MESSAGES 

The following error messages are typed at the terminal 
as soon as they are detected. If the error is only in 
syntax, the system is immediately ready to accept a new 
query, otherwise execution is terminated. Errors have been 


grouped with one message for each group. Messages are selt 


explanatory. 
ERROR NO POSSIBLE REASON 
3. >: Attempt to parse a blank query. 
OL : Diack cetyepmtiela. One of the characters 


DV ,¢ of blank was expected after A- 


oy 





102 


iOS 


104 


201 


Z02 


301 


302 


208: 


904 


S105 


506 


= O1%, 


308 


Sal 


Invade ivyremereld. "One of “the Charcacthens 
QO,E was expected after S. 

[Invamwed TYPE iveld> One ot the characters 
P,V was expected after NO. 

Invalid TYPE field. No command starts 

with the given letter. 

Missing character or somewhere in the 

query there is no space delimiter. 

There 1S no space delimiter. 

Invalid SELECT field. SELECT field is 
missing or there iS no space delimiter 
between TYPE and SELECT fields. 

Invalid CONDITION field. the word FOR is 
missing (the string OR was expected after 
F), or invalid PRINT field. 

Missing left parenthesis in CONDITION field. 
Incempletce condition “field or missing 
space delimiter. 

Missing right parenthesis in condition field. 
Invalid OR logical operator. Character R 
R was expected after OQ. 

Invalid AND logical operator. The stran¢ 
ND was expected after A. 
inValrdelogicalsoperater.  @mly OR” and 
AND are accepted. 

Invalid operand for statuS. The charactes 


C or S was expected after B. 
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5) 


3 


514 


=a 


Skis 


Sey 


eS 


Se, 


570, 


601 


602 


603 


604 


605 


Invalid operand for status. The 
L was expected after L. 

Invalid operand for status. The 
L was expected after U. 
invealtadsoseband £or Status. The 
Q or R was expected after E. 
Invalid operand for status. The 
Vo0 5) ao en 5 a expected before 


Invalid operand for status. The 


character 


character 


character 


character 


on 


character 


Ver S pi, wi, -fmes, | Was Gsspecteds betenre Rk. 


Invalid operand for status. The character 
P was expected before P. 

invalad Gperand for status. The character 
I, A, C or D was expected before E. 


Missing space delimiter after status 


operand. 


Non recognizable operand for status. 


inialoecderLirst operand for relational 


@perator in Cconditlenmmrr1eld.: 


Missing space delimiter in a simple 


conditional phrase. 


Invalid relational operator. The character 


T or E was expected after G. 


Invalid relational operator. “The Character 


T or E was expected after L. 


Invalid relational operator. The character 


Q was expected after E. 


oo 








606 > Invalid relational operator. The character 
E was expected after N. 
607 : Invalid operand in bound conditional 
phrase. The character N was expected 
after string MI. 
608 3 invalid operand in bound conditional 
phrase. The string AX waS expected after M. 
609 : Unrecognizable relational operator in 
Simple conditional phrase. 
701 Teinvalela prince t1eld, or missing word FOR 
Ln econemeron field: 
1001 ; Error in input data. Unrecognizable 
status code. 
1002 : Error in input data. Data encountered 
has less than the expected number of rows 


and columns. 


Pees METATIONS —- EXTENSIONS 

As mentioned in the introduction, ATHENA is a direct 
expansion of the PERUSE System. It includes all the features 
Of PERUSE, except the weighted average command and has the 
following differences and extensions: 

1. ATHENA supports two distinct data structures, each 
different from that of PERUSE. This was neceSSary in order 
to support efficient access to individual records or group 
of records. The SUPERSPARSE data structure 1S unique to 


ATHENA. 
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2. ATHENA accepts as input a simple file which can be 
easily obtained from the solution file of any linear 
programming package on tape, disk or cards. 

3. ATHENA supports the commands SET, ACTIVATE, 
DEACTIVATE and COUNT, in addition to the commands of 
PERUSE, allowing the user to construct logical subsets of 
the solution and efficiently access these subsets as inde- 
pendent files with very small access time. 

4. ATHENA supports compound conditional phrases for 
extraction of more specific information and the bound 
conditional phrase for maximum and minimum values. 

59. ATHENA uses object time variable format allowing 
better appearance of output and uses the words NONE and 
INFINITY instead of the number 0.7273E76 for better 
readability. 

6. ATHENA accepts reduced problems from PREP [2] and 
can be used to pass the PREP status file with the solution 
file of any optimization system to permit recovery of the 
Original problem solution. 

ATHENA has been designed to handle solution files with 
up to 30,000 records. The actual limit 1S imposed by the 
number of real number values that must be stored explicitly. 
This number cannot presently exceed 65536 since this is the 
largest integer pointer value which can be stored by ATHENA 
in a 2-byte halfword. Experience has shown that the average 


number of stored values for each record, excluding penalties 
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meee i7/). Adding to this another 0.5 per record for 
penalties, the problem size limit may be as large as 30,000 
rows and columns. | 

A rough estimation of the space needed for the packed 
file in 4-byte words can be obtained by multiplying the sum 
of rows and columns of the solution file by 5 for the SPARSE 
data structure and by 4 for SUPERSPARSE. Before using 
ATHENA, adjust the size of the SOLFIL array in common block 
SOLPAC to this number. To avoid passing problems with common 
areas under the IBM 360/67 Operating System, use an array 
size that 1s an exact multiple of 4096 larger than the 
number calculated above. Also make corresponding adjustments 
to the DEFINE FILE statement of the main program. Other 
computers will require analogous modifications to this array 
size. 

ATHENA has been developed in modular form and can be 
easily changed or extended to support future needs. Commands 
which can be easily implemented include the weighted 
average, the sort of output, or further calculations needed 
for the analysis of the solution file. ATHENA can also be 
used aS part of an integrated system for sensitivity analysis 


Of optimization problems. 


ie ofoltEM INTERFACE 
1. Preview 
Linear Programming packages give differing forms of 


Output so that itis difficult for a system to be interfaced 
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with all of these solution formats. ATHENA accepts as 
input a solution file in a 'standard' form which can be 
easily obtained from any other solution file form. 

ATHENA is best utilized in an interactive system 
such as CP/CMS, although it can also be used in batch 
processing. On the other hand, most L.P. packages run only 
in batch processing. Moreover, in some systems there is 
no integration of interactive and batch processing (as is 
Surrently the case at the Naval Postgraduate School Computer 
Center). In these cases, the solution files may be transferred 
manually from one system to the other using magnetic tapes 
eee cards. 

A simple input file has been designed which can 
be punched in cards or entered on tape, disk, or other 
storage media. 

ee enputc File 

Picowinpue rit le Consists of records with the following 
Structure: 

a. The first record always contains the number of 
ROWS, the number of COLUMNS, and in position 51 the character 
'1l' if each ROW record contains PENALTIES or '0O' otherwise. 


The FORMAT of the first record is 
GneresO x, 15440 Xl). 


b. Each subsequent record contains explicitly all 


the information associated with each ROW and COLUMN, with 
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the following format: 


NAME Format — 2A4 (left justified) 


STATUS ROrmat A2 


mec or S, L, U, and D numeric field values with Format 
5E14.5 or 5F14.5. (The meanings of each of these fields is 
described in the previous section.) If the solution file 
includes PENALTIES, then two records will be associated 
with each ROW. The first will be exactly that described 
above, and the second will have the FORMAT (E14.5,16X,E14.5) 
eee, and W, the UPPER and LOWER PENALTIES. In all cases, 
INFINITE values will be represented explicitly by the number 
+0.1E76. The total number of records must agree with the 
sum of ROWS plus COLUMNS, with the records of ROWS preceeding 
those of COLUMNS; otherwise an Input error will occur. The 
file 1s read in and packed one record at a time. 
3. Packed File 
a. Packed File as Input 

If the input file is already packed from a 
previous use of the system, it will be read in unformated 
binary format. The system will provide the user information 
for memory requirements before reading the file. The packed 
file may be on a tape or disk but cannot be on cards. The 
system will ask the user at the beginning of a session for 
the number of the file. 

b. Packed File as Output 
If an unpacked file is used as input, the system 


will ask for the file number where a packed file is to be 
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written. Of course, a DEFINE FILE FORTRAN statement must 


be included right after the declarations of the main program. 


meee SAMPLE OF SYSTEM USE 

The procedure follows for use of ATHENA at the Naval 
Postgraduate School Computer Center with the IBM 360/67. 
The solution file here is produced by the MPS/360 package 
and ATHENA is used under CP/CMS. Similar procedures can be 
followed for any other installation. 

ie Obtaining the Unpacked Solution 

a. Submit the problem to be solved uSing the usual 

Control Cards required for the MPS/360 package inserting 


before the Control Card: 
//MPS2.SYSIN DD * 
eave £Ollowing cards: 


J//MPS2.SYSPRINT DD DSN=Sxxxx.nnnnnn, 
i ont r=3330, 701l=SER=DISK04, 
// SPACE=(CYL, (1,1)),DISP=(NEW,KEEP), 


// DCB=(RECFM=UA, BLKSIZE=133) 


With these cards the output of MPS will go to the DISK 


instead of the printer. 


XXXX 1s the user's number and 


nMnnnnn is the file name on the disk. 
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Pe etemcOiMtaon £ile 15 t00 big or disk epace 
is not available use tape or tapes to store the output. 

b. Use the program REWRITE (page 98) to transform 
the MPS/360 tape or cards to the format required for the 
ATHENA unpacked Input file. 

c. Now the unpacked file for ATHENA is available 
and can be used to analyze the solution. 

ae Using the System Under CP/CMS 

ATHENA in CP/CMS TEXT form requires about 67K bytes. 
The space needed for the packed file depends on the number 
of records, the method used for packing and the density of 
the original file. 200K bytes would be sufficient to hold a 
packed file with up to 13,000 rows and columns. After 
Sufficient space has been secured, the following procedure 
may be applied: 

a. Ask OPERATOR to connect the tape with the Input 
file created by the REWRITE to the private disk as device 
181. As soon as the tape is connected, the message 'DEVICE 
181 ATTACHED' will be printed at the terminal. 

b. Before using the tape, type ALWAYS under CMS 
the command ' TAPE SKIP 1 '. This command will position 
the tape at the first record of the file. This command is 
required because the tape created by IBM 0.S./360. 

c. Type $ ATHENA 

ATHENA will ask for information about the file 
identifiers for input-output, whether the file is packed 


and method of packing and will give the size of the packed 
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file. For input file ordinal use any number between 01 
ea 99 excluding the numbers 03-06. For output file 
ordinal give the number 03 or 04. These are the numbers 
used by the DEFINE FILE FORTRAN statement and they can be 
changed. The system will be ready to accept QUERIES as 
soon as the prompt phrase 'INPUT A COMMAND' is typed by 
ATHENA at the terminal. The packed file can be saved on 
a tape uSing the 'TAPE DUMP' command under CMS. 

In the next few pages a demonstration of using 
ATHENA with a solution file of 766 rows and 10921 columns 
is given. This problem is a mixed integer optimization 
model with 963 binary variables for medium term capital 
budgeting of the Naval Air Test Center [6]. For this 
problem, a query may require as much as one and a half 
Minutes of clock time if the interactive system on the IBM 
360/67 1S under heavy yes and the query is difficult to 
answer. However, most queries are answered almost immediately. 
Response time is especially good when the user makes use of 
ACTIVATE, MASK and SET features to qualify necessary 


searching. 
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$ ATHENA 
Pea) LS THe FILE NO OF THE SOLUTION FILE ” (FORMAT I2) 
4 
153 THE FILE ALREADY FACKED ? ENTER ‘YES OR NO 2? 
iS 
menoey REQGUITRMENTS FOR THe FACKED FILE sAZ 7A = 4-EVYTE 
WORTIS 


meee Wave SUPP ICTENT MEMORY SPACE ENTER YES 


DQTHERWISE ENTER NO MAKE ADJUSTMENTS AT COMMON ANT TRY 
AGA TN 
ris 

TNFUT & COMMANT 

MOM KCK HOH HCHO TK MOH SHAK CHCA HO HK MO HEHE CAH CHO HHA HK MK OK 
Poor oe fhe CHeERACTER * AS THE FIRST COLUMN OF THE QUERY « 
THE QUERY IS TGNORED ANU THIS [35 A CONVENIENT WAY TO x 
Peoear COMMENTS I TRE QUIFUT FROM A TERMINAL. 
COMMENTS IN THE QUTFUT FROM THE OFFLINE FRINTER ARE 
PeoeRTeEl USING THe LETTER H INSTEALR OF k. 
TO AVOLD THE FROMFTNESS FHRASE “INPUT A COMMAND’ AT 
we OUTPUT VSe THE COMMAND *“NORPROMFT’. 
Recess: GU THe REQUIREMENTS FOR THE APPEARANCE QF THE 
Maeots (HE FRINT G@EPTIGMN FIELD OF THE UISFLAY COMMAND x 
1K feewoseh ONLY WIThH AT MOST TWO OFTIONS. ¥ 
MOSHE HONG KOK CK CHOSE K HOM HOH MH HK HK CK Me HR Me MEE GK HK KK OK 
4 
Rd 
K HERE IS A GEMONSTRATION OF USING THE SYSTEM ATHENA. 
x 
Me 
¥ HOW MANY ROWS HAS THE FILE ? 
K 
COUNT ROWS x 


Soot Se SS NK Se 
~~ o> OT OX*« 


ray 


coc 


765 ROWS QR COLUMNS WITH MASK ¢ KKK KKK 
Seep kSsrY THE CONDITIONS 
x 
eK HOW MANY COLUMNS °? 
BY 
COUNT COLUMNS <x 


LOV2SL ROWS OR COLUMNS WITH MASK ¢ k:koCKKKEK 
Petar yY THE CONUITTIONS 


kK HOW MANY OF THEM ARE BASIC 7 
x 


te fee te FORCST &S) 
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504 ROWS OR COLUMNS WITH MASK > tcKxOKaCKKx 
SATISFY THE CONDITIONS -- 
k HOW MANY QF THEM ARE EQUAL ZERO ? 
k 

CC * FORC X EQ oO. AND ST BS) 


ga ROWS OR COLUMNS WITH MASK 2 KKKKKKKK 
Ser lsreY THE CONDITIONS 


— 


PeiNG THE COMMAND “ACTIVATE’ ONLY THE ACTIVATET 
SeeukUS ARE SEARCHED TO ANSWER THE QUESTION. 
mers tom GUCD WAY TO AVYOTT SEARCH OF THE WHOLE 
ae. 


* KR KK 


we 


i 
ACTIVATE COLUMNS XO 


Es Pio s weet TVA Tr i 
poral RECORDS ACTIVE ; os 


*K HOW MANY ROWS NOW 


O ROWS OR COLUMNS WITH MASK AK RHORKE 
Siri ioryY THE CONDITIONS 


fence Fike NO ROWS SINCE ONLY COLUMNS ACTIVATED. 
Dow Meavy OF THEM ARE BASIC * 


ce eo 


See * FORCST 8S) 


Ie ROWS OR COLUMNS WITH MASI 3 HORROXKOKHKX 

Sai lsreyY THE CONDITIONS 
x 
x Pees lar THEM 


es, 


7 


MeerlAy C xX FORCST BS) XxX L. 
Seer WELOWING ROWS OR COLUMNS SATISFY CONDITIONS 


NUMBER  .NAME.» AT «s ACTIVITY... eee yal sles ee -« 
B03 XO4d eve Oy reo O9,O008D0 


1 ROWS OR COLUMNS WITH MASN 3 KKKXCKRKKK 
SATISFY THE CONUITIQNS 
K 
XK WHAT IS THE AVERAGE OF NONZERO ACTIVITIES ? 
K 
AVERAGE C * FORCX NE O) 
THE FOLLOWING ROWS OR COLUMNS SATISFY CONDITIONS 
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NUMBER NAME... AT «--ACTIVITY... pommeeT COST. . 
Mero OK AVERAGES ¢ . ee Oe oad oS 


71 ROWS QR COLUMNS WITH MASI KKKKKKKK 
Speeter yY THE CONDITIONS 


x 


Ree sUM OF ACTIVITIES AND INFUT COST 7 


ras 


meu © * FOR ¢ X NE O 3 
eee ULLOWING ROWS OR COLUMNS SATISFY CONDITIONS 
Borer NAME... AT »..-ACTIVITY... ee eae cs) fey 


SUMS OR AVERAGES ¢ Poe eGo? J%.363882 


71. ROWS OR COLUMNS WITH MASK 3 OKA 


Mees Y THE CONDITIONS 


PaeeoomMa NEW COLUMNS AT THE ACTIVE FILE 


> 3 3 


mee 86YMOL 


_—, 


Tea RECORDS ACTIVATED 
Retr RECORUIS ACTIVE ¢ 243 


Biot F 


Peer UR (ST RS) 


7 ROWS OR COLUMNS WITH MASK §  ACKKOK ARK 
eel iskyY THE CONDITIONS 


af 
*K feeor bar BASIC AND ZERO ACTIVITY 
x 


Ome 6 KUL FLOOR UC UXUC EQ OO AND STATUS BS) xX C 
Peer ULLOWING ROWS OR COLUMNS SATISFY CONDITIONS 


NUMBER »>NAME * + AT > @ Se) Ey a *» INFUT Best ¢ > 
moos  YMOLOdéL Bs 0.00000 SO) iO Or are 


1 ROWS OR COLUMNS WITH MASK 3 AKA 
eeailsry THE CONDITIONS 
x 
K pewetiining INPOUP COSTS ARE ZERO * 
. 
mee x FOR C EQ 0 ) 
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% ROWS OR COLUMNS WITH MASK 2 KRKMKKH XK 
GATISFY THE CONDITIONS 


x 
iK HOW MANY AT LOWER LIMIT * 
XK 

Meee ® FORCSTATUS LID 


171 BOWS OR COLUMNS WITH MASK & KKKKRKKK 
Semeoty THE CONDITIONS 
x 
x Sethi: SOME RECORDS FROM THE AUTIVE FILE 
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APPENDIX 


The whole program consists of 3 main parts: The input 
part, the parser-code generation part and the interpreter 
part (see Figure 1). The main program calls these 3 parts 
according to the instructions from the user. Mcre detailed 
information for the main subroutines of each part are given 
below. 

im 6CLNpDUt Part. 

There are 2 almost identical subroutines named 
READF and READD. READF is uSed to pack the solution file 
in SPARSE data structure and READD in SUPERSPARSE. 
a. Input Parameters 
(1) IFLR: The file number of the solution 
file. This can be any number between 1-99 excluding the 
numbers 3, 4, 5, and 6. 
(2) IFLW: The file number where the packed 
file will be written. Currently this number can be 3 or 
4 since the DEFINE FILE statement in the main program defines 
the files 3 and 4 only. 
b. Output Parameters 
(1) NROWS, NCOL; The number of rows and columns 
of the solution file respectively. 
(eer kR: ff chis value is different from zero 


an error has been encountered in the input data. 
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The solution file iS read in standard format, 
as explained in user's manual, one record at a time, filling 
the array SOLFIL in the common area SOLPAC with the packed 
file. The name of the record is read directly into SOLFIL 
and the associated values into the array DATA. The 2 
character status code 1s read into the INTEGER*2 variable 
STAT. The INTEGER*2 array CSTAT has all possible values for 
Status, and their location in this tesa Minus one, gives 
the corresponding value of the first bit group in the 16 bit 
map aS explained in the data structure section. 

The variable INDEX always points to the next 
avallable location for storing the next nonzero, noninfinite 
number. The value of INDEX at the beginning of the processing 
of each record is stored in the 2-byte pointer after the 
name. Each number in the array DATA is examined and if it 
is eligible for storage it is stored at the location 
SOLFIL(INDEX) and INDEX is increased by one. At this point, 
for the subroutine READD, the end of the array SOLFIL, where 
the distinct real number values are stored, is searched and 
if the value sought is found the location of this value, 
counting from the end of the array is stored at the location 
SOLFIL(INDEX); otherwise, a new entry is made and the location 
Seeenis entry is stored at the location SOLFIL(INDEX). Notice 
mae SOLPIL(INDEX) is a 4—byte location for READF but 2-byte 
for READD as explained earlier. After the last record has 


been packed, READD shifts the distinct real number values 
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down to the end of 2-byte indirect pointers. When the 
System is to be used under batch processing the CALL TAPSET 
statement from READF and READD must be removed and the 
input file must be described using the usual DD statement. 
2. Parser = )Code Generator 
There are several subroutines in this part of the 
System. There 1S one subroutine for each field of the 
query, the control subroutine PARSE and various others 
to perform specific tasks like getting the numerical value 
of a string of decimal characters,etc. The subroutine PARSE 
is described below. The numerical codes used for each case 
have been inSerted as comments at the beginning of each 
Subroutine where they are used. 
Gee Lnpuc Pagamecers 
(1) QRY: INTEGER*2 array with the query to be 
parsed in A format. 
b. Output Parameters 
(LjeeGebE: INT ECER*2,array with the numerical 


code for query as follows: 


ORY (1) wey PE Piledd 
RY (2) Peok LECT field 
ORY (3) > MASK field 
eey (4) - ORY(10) : CONDITION field 
Sewell) - ORY(18) : PRINT field 


(2) MASK: INTEGER*2 array with the mask from 


the query, one character per location. 
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(3) ARG2, ARG4 : numerical values from the 
menattion field, if any. 

(4) Ph code numMeer £LOr SYMntax error. If 
TER=0 the query has been parsed correctly. 

Subroutine PARSE calls sequentially the corres-~ 
ponding subroutines to parse each field of the query. Before 
calling it positions the pointer ICOL1 at the location of 
the QRY where the field starts using the subroutine GETNCH. 

Pepe meen pLeeter GPAre 
There are several subroutines called by the main 
Subroutine INTERP of this part. 
ar oneuc Parmacters 

CE eCOUE,. MASK; RG2z, ARG4, NROWS, NCOL, IER 
as described previously. 

(ee ne ee RecOrd number to start 
Searching, to stop searching and searching step respectively. 
Default values are 1, (NROWS + NCOL), l. 

(3) IFLG : Flag to determine the data structure 
type of the packed file. If IFLG=0, data structure SPARSE 
Otherwise SUPERSPARSE. 

(4) IACT : Passed through the INIT common area 
and, if 0, the searching is done as determined by Kl, K2, 
K3, otherwise there is an active file created by the user 
through the ACTIVATE command andonly the members of this 
file are searched. 

The logic of this subroutine is shown in the 


flowchart on the next two pages. 
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